Magnetic penetration depth in the presence of a spin-density wave in multiband 

superconductors at zero temperature 
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We present a theoretical description of the London penetration depth of a multi-band supercon- 
ductor in the case when both superconducting and spin-density wave orders coexist. We focus on 
clean systems and zero temperature to emphasize the effect of the two competing orders. Our cal- 
culation shows that the supefluid density closely follows the evolution of the superconducting order 
parameter as doping is increased, saturating to a BCS value in the pure superconducting state. 
Furthermore, we predict a strong anisotropic in-pane penetration depth induced by the spin-density 
wave order. 

PACS numbers: 74.20.-z, 74.70.Xa, 74.25.N-, 74.25.Dw 



o 

o 
i 

5h 
Oh- 
3 ■ 

! 



o 
o 



> 

in 

in 
o 

(N 



X 



I. INTRODUCTION 

Most iron-based pnictide superconductors exhibit a 
phase diagram with the spin-density wave (SDW) and 
superconducting (SC) phases neighboring each otheri^ 
More detailed studies of the boundary between the 
two phases indicated that the SDW and s" 1 SC order 
parameters^ may coexist on a microscopic level^— due 
to the two-dimensional multi-orbital character of elec- 
tronic band structure of these materials^— Measure- 
ments of nuclear magnetic resonance^ and neutron scat- 
tering experiments 7 - on electron-doped BaFe2As2 com- 
pounds support the presence of the coexistence phase 
between a pure superconducting and a pure SDW phase. 

The competition between the SDW and SC orders 
tends to suppress one order as the other order devel- 
ops. Analysis of the relative weight of these two antag- 
onistic order parameters reveals the microscopic nature 
of the coexistence region of the phase diagram, which 
can be assessed, for example, from thermodynamic char- 
acteristics of materials. Indeed, as doping decreases 
and the material moves towards a stronger SDW phase, 
the specific heat jump at the onset of superconductiv- 
ity decreases ii2r— . Another thermodynamic quantity 
that shows a suppression of the SC order by the SDW 
ordering is the magnetic penetration depth. Its measure- 
ments in iron-based superconductors is focus of many 
experiments (see references provided in Ref. 0) , and re- 
cently were reported in the coexistence phase of Co- 
doped BaFe 2 ^ a: Co a ;As2 compounds! 16 ' 17 

The London penetration depth A manifests itself as a 
length scale that determines the penetration of a static 
magnetic field within the bulk of the superconductor. 
The magnetic field inside the superconductor is given by, 
see Fig. Q^a): 



H(y) = H ext exp (-y/X) 



(1) 



At zero temperature, A is defined by the supcrfluid den- 
sity n s in the two-fluid model^ 



1 



n s q 



(2) 



where q is the charge, and m* is the mass of the charge 
carriers in the supcrfluid phase. The BCS theory pro- 
vides the same result, with the caveat that the charge 
carriers are Cooper pairs that carry charge q = 2e, and 
have a mass m* = 2m. This result for A does not de- 
pend on the actual strength of the SC order parameter 
at z:ero temperature when n s is the total electron den- 
sity. Because electron doping is not expected to drasti- 
cally change the Fermi surface of BaFe2- a; Co a ;As2 com- 
pounds, the observed increas e) 16 ' 17 of the magnetic pen- 
etration depth at zero temperature is likely due to the 
competition of superconductivity and SDW. 

In this paper we present a microscopic theory of the 
magnetic penetration depth in the coexistence phase be- 
tween SDW, and s H SC statesJ^ We utilize a simpli- 
fied model for iron-based pnictides as two-dimensional 
multiband metal. We focus on the zero-temperature limit 
when the magnetic penetration depth can be interpreted 
as the density of quasiparticles participating in formation 
of the SC order, while finite temperatures lead to the re- 
duction of the superfluid density, and increase in the Lon- 
don penetration depths We also disregard the effects 
of impurities on the magnetic penetration depth. There- 
fore, we isolate the effect of the competition of SDW with 
the SC order parameter for the available electronic states 
from other known effects on the magnetic penetration 
depth. 

Within the above assumptions, we obtained the mag- 
netic penetration depth A in the coexistence phase. We 
show that A is determined solely by the density of quasi- 
particles related to the SC order parameter. Thus, as 
the doping moves the system from a pure SC phase into 
the coexistence phase, an SDW order develops and takes 
quasiparticles away from the SC condensate. As a conse- 
quence, A monotonically increases and becomes infinite 
at a transition to a pure SDW phase. This behavior is 
consistent with the results of measurements of the mag- 
netic penetration depth in the coexistence phasei 16 ' 17 

Our model also predicts an anisotropy of the magnetic 
penetration depth. Earlier theoretical model a 10 ' 11 pre- 
dicted that the coexistence between homogeneous SC and 
commensurate SDW orders requires a weak anisotropy of 
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FIG. 1. (Color online) a) Exponential decay of magnetic field, 
(a, fo, c) are the principal axes of the sample. The 1/e decay 
length is determined by the Q a principal value of the Q-tensor. 
7r labels the SDW stripe order vector; b) electronic structure 
of the multi-band model in the unfolded Brillouin zone. The 
hole FS is centered on the V point, with SC order parame- 
ter Ah, and the electron FS pockets are centered at (±7r, 0), 
(0, ±7r), with SC order parameter A e . The magnetic order 
m with momentum 7r = (tt, 0) hybridizes hole and electron 
states separated by tt, but leaves FSs at (0, ±7r) intact. 



an electron or hole Fermi surfaces (FSs). We demonstrate 
that even this weak anisotropy is sufficient to produce 
a significant in-plane anisotropy of the magnetic pene- 
tration depth. At its peak value, the anisotropic part 
is about 30% of the value of the isotropic part. Unlike 
the isotropic part, this anisotropic part manifests itself 
within the coexistence region and vanishes in a pure SC 
phase. Measurements of the London penetration depth 
without taking into account the relative orientation of 
the crystal interface may lead to a large scattering of 
measured A in samples with identical doping, observed 
in Ref. H 

The outline of the rest of the text is as follows: Section 

II is a summary of the main results, together with a short 
description of the model within which we work. Section 

III contains the technical details of our calculations of the 
magnetic penetration depth in the coexistence phase. We 
present our conclusions in Section IV. 



II. MAIN RESULTS 

We consider a multi-band modeli^iii of iron-based 
pnictide superconductors as an itinerant metal contain- 
ing two electron, and one hole bands, see Fig.QJb). The 
FS for the hole band (h) is centered at the T-point {0, 0} 
of the unfolded Brillouin zone and the FSs of the elec- 
tron bands (ei 2) are centered at the M-points, defined 
by vectors 7Tx = {±7r,0} and 7T2 = {0,±7r}. At perfect 
nesting, the band dispersions satisfy £ e (fc) = —£h{k), 
where k is measured from the center of the correspond- 
ing band, and the two FSs are circular with equal sizes Uf 
(kp <C tt/2). A deviation from a perfect nesting may be 
tuned by doping or pressure and is parametrized by the 
small detuning parameter 5(k) = (£ e (fc) + £h(fc))/2. An- 
other quantity that enters in later calculations is £(fc) = 



(£ e (fc) — £h(fc))/2. At detuning from perfect nesting, the 
hole pocket remains circular, whereas the electron pock- 
ets may change their size, as well as become elliptical. 
Both of theses effects can be accounted for by the ex- 
pression 

Sf, = S + 8 2 cos(2</>), k = {cos(0), sin(</>)} (3) 

with two band-structure parameters 5o,2- We note that 
we need to consider ^0,2 ^ T s _q <C vf kp, where vp is 
the Fermi velocity of the fermions at the FS. We treat the 
two-body interactions in the SC and SDW channels in the 
self-consistent approximation by introducing SC, A, and 
SDW, m order parameters. We assume that the super- 
conducting state has an extended s H -wave symmetry^ 
and is characterized by an equal in magnitude, but oppo- 
site in sign isotropic SC order parameter on the FS, al- 
though a more general case, that also includes anisotropy 
of SC order parameter— >2i deserves a separate consider- 
ation. The SDW order parameter is characterized by 
the nesting vector tti^ that establishes an SDW stripe 
structure and breaks the C4 crystal symmetry^ The 
self-consistency equations for A and m are obtained from 
those derived in Ref. [ll| for finite temperatures by taking 
the limit of T — > 0, as outlined in the end of Section III. 
The equations are^ 



In- 



dcu Re 



(E + iS^/E 



(4a) 



In 



[ c,0 



T s ,o Aq 



duj Re 



y/(E + i6tf + H 



where E = + |A| 2 , and (...) stands for averaging 
over FS. Within this model, the SC transition tempera- 
ture from a normal metal to a pure SC phase does not 
depend on detuning and is labeled by T Cj o, while the tran- 
sition temperature from a normal metal to a pure SDW 
phase decreases as the detuning from perfect nesting in- 
creases. To characterize the interaction strength in the 
SDW channel, we introduce the transition temperature 
to a pure SDW phase at perfect nesting as T Sj o. Param- 
eter Ao in Eqs. ((3]) is the SC order parameter at zero 
temperature for a pure SC phase, Ao = 7re -7E T Cj o, with 
je ~ 0.577 being the Euler constant. 

The expressions for the response function Qik that con- 
nects the supercurrent density and the vector potential 
as ji = — QikAk are derived in the next section. The 
tensor Qik can be written as a combination of isotropic, 
<5 ISO , and anisotropic, Q an , components that are given by 



Q" 



Q A 2 f d£d<i> 
2tt 



1 

cos 2(f) 



£ 2 Y(0 + 4m 2 Z{£) 



e 



(5) 



at the zero temperature limit. Here 
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2E+E-(E+ +E-) 



e±(z, 0) - ^(v^+e ± s h r + a 2 

and Qo is given by the expression 
„ 2e 2 t;|JV> 

yo = • 



(6) 



(7) 



(8) 



in terms of the Fermi velocity, vp, and the density of 
states at the Fermi energy, Np, and does not depend on 
configuration of states far away from the Fermi surface. 

The electromagnetic response tensor Qife has the prin- 
cipal axes that coincide with the crystallographic axes in 
directions tti^. If the crystal surface is perpendicular to 
the iron plane and is along one of these directions, the 
magnetic penetration depth have simple expressions 



\Q a 



4?r A 2 

max 



Q iso -\Q an \. (9) 



In view of these simple relations between the magnetic 
penetration depth and the response tensor Q, below we 
discuss the properties of the electromagnetic response 
tensor. 

In Fig. [21 we present a case with parameters 
T sfi /T cfi = 2.0, and 8 2 = 0.4(2ttT c , ). Panel (a) dis- 
plays the phase diagram^ as the isotropic detuning pa- 
rameter Sq is varied. The filled area represents the re- 
gion of coexistence between SC and SDW orders. As 
the temperature is increased, the range of values for So 
where coexistence is possible shrinks, degenerating into 
a tetra-critical point. We note that within this model, 
detuning plays a detrimental role only for the SDW or- 
der parameter. This is clear from the phase diagram, as 
we can see that the Neel temperature decreases. Panel 
(b) presents the evolution of the zero temperature order 
parameters A/Ao, and ra/mo, normalized to their val- 
ues in pure SC, Aq, or SDW, tuq, phases, respectively; 
note that Ao/T c .o = mo/T s< o = 7rexp(— ^ E )^L As the de- 
tuning increases, the SC order parameter monotonically 
increases from zero at the lower edge of the coexistence 
region to the full BCS value at the right edge of the co- 
existence region. The SDW order parameter follows an 
inverse trend. Using the values of the order parameters, 
the zero-temperature value for the response tensor is cal- 
culated according to Eq. ([5J. We note two limiting cases. 
First, when A = 0, the integrand for both the isotropic 
and anisotropic component is a finite number, and, be- 
cause of the prefactor A 2 , both Q ISO and Q an vanish at 
A = 0. Second, when m = 0, the integrand decomposes 
into a sum for two independent bands ±£ + 5? . By a 
change of variable, the integrand can be made isotropic. 
Then, the integral for the isotropic component reduces to 
1/A 2 for each band, and we obtain the value Qq, Eq. ©. 
The average over the angles for the anisotropic compo- 
nent in Eq. (J5J) is zero due to the cos(2</>) term. These 
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FIG. 2. (Color online) a) The phase diagram of a two band 
metal with the interactions in SDW and SC channels as a 
function of the detuning of the FSs; b) The zero-temperature 
SC and SDW order parameters as a function of detuning So 
of the FSs; c) The isotropic (solid) and anisotropic (dashed) 
components of the Q-tensor as a function of detuning of the 
FSs. T s ,o/T c , = 2.0, and 5 2 = 0A(2ttT c , ). 



two limiting cases determine the behavior of the plots 
in panel (c). The isotropic component rises monotoni- 
cally from zero to its BCS value, whereas the anisotropic 
component goes through a maximum in the coexistence 
region, and becomes zero at both edges of the coexistence 
phase. 

We present the minimal and maximal values for Q 
in Figs. [3] and 0] for two cases: a) T S:0 /T C:0 = 2.0, 
and 6 2 /(2irT c . Q ) = 0.4, and b) T S)0 /T Cj ' =' 5.0, and 
82 1 '(27rT C: o) = 1.0. The difference between the plots is 
that in Fig. [3j Q is presented as a function the dimension- 
less isotropic detuning parameter So / (2nT C fi) , whereas in 
Fig. 01 Q is plotted versus the transition temperature, T c , 
from a pure SDW phase to a co-existence phase, obtained 
for the same FSs detuning parameters <5o,2- Q is normal- 
ized in units of Qo, Eq. ([5]), and T c in units of T c , which, 
in our model plays the role of the highest critical tern- 
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FIG. 3. (Color online) The maximal and minimal values 
(solid) for the Q parameter, and the value for no detuning, 
Eq. (|10l) (dashed) for different values of the ratio T St o/T c ,o, a) 
T s fi/T c fi = 2.0. The coefficient for the anisotropic component 
of the detuning is S 2 = 0.4(2ttT c ,o); b) T sfi /T cfi = 5.0, and 
S 2 = 1.0(27tT c , ). 

perature reachable at optimal doping. Figure U attempts 
to eliminate model-dependent detuning parameters #0,2 j 
and express the dependence in terms of quantities that 
are available from direct measurements. In all cases, the 
presented dependence is a monotonic function. However, 
we observe that the anisotropy is relatively strong in the 
coexistence phase, which is manifested by a large spread 
of the bands formed by these plots. As the strength of 
magnetic coupling relative to SC coupling is increased, 
represented by the increase in the ratio T s q/T c q from 
2.0 to 5.0, the spread of these bands increases. This is in 
line with the notion that the anisotropy is caused by the 
SDW order, which scales as T 8) o- 

The plots on Fig. @]bccome flat at the right edge of the 
coexistence region. This flat region is understandable if 
one considers the phase diagram in Fig. HJa). Namely, 
the right end of the coexistence region shifts towards 
higher values of the detuning parameter 6q as temper- 
ature increases. Thus, the tetracritical point takes place 
at higher detuning 6q than the value of <5o for transition 
from coexistence SDW+SC phase to a pure SC phase 
at zero temperature. In this window of Sq, the response 
tensor is fixed, Q = while T c still monotonically ap- 
proaches T c fi, This gives rise to the flat part of the curve. 
Notice that this flat region is significantly reduced as ra- 
tio T Si q/T Ci q increases, since the above mentioned window 
of values for Sq shrinks as welliii 



Finally, we note that if we disregard the detuning pa- 
rameter in expressions for the response functions, we ar- 
rive to the following simple expression^ 



The co-existence between the SDW and SC order param- 
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FIG. 4. (Color online) A parametric plot of Qmax/Qo 
(dashed) and Q m i n /Qo (solid) vs. T c /T Ct o, as the detuning 
of the FSs is varied for different values of the ratio T St o/T Ci o 
a) Ts.o/Tc.o = 2.0, and coefficient for the anisotropic compo- 
nent of the detuning S 2 = 0.4(27rT c , ); b) T s , a /T c . = 5.0, and 
h = 1.0(27rT c>0 ). 



eters requires a finite detuning with non-zero values of i5o 
and 62 , and the detuning between the FSs has to be taken 
consistently both for evaluation of order parameters A 
and m, as well as for calculations of the response func- 
tion Q. Dashed lines in Fig. [3] represent the result for Q 
according to Eq. (fT0|) for order parameters evaluated for 
zero temperature limit. This estimate for the response 
function is smaller than even the minimal value of the 
anisotropic response function obtained from Eq. (JS|). 



III. CALCULATION 

Here we work within a model developed in Ref. [ill For 
definitcness, we consider the case of coexistence between 
an -wave SC order, and a commensurate SDW order. 
The system is defined through the mean-field Hamilto- 
nian: 
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ft = | Efc Q( 3 *L (ftfcW *fe|8 + Ene{/»,e},fc £ri(k), 







m *( d ^ a B 





-A* h (ia 2 )af3 


-£h(-k)5 a p 












£e(k)<W 







-m( ( T 3 r ) aP 


-A*(i<T 2 ) a ^ 





with ^ ka = (c\ ha , c h - ka , c\ ka , c e - ka ) and ^ ka being the 
conjugated column. c eka and c^q, are electron and hole 
band annihilation operators, respectively. The two diag- 
onal blocks of the matrix (ftfc) Qj a correspond to a purely 
SC system with Ah and A e developing on two different 
bands, and two off-diagonal blocks contain SDW field 
m = m z that couples fcrmions between the two bands. 
In the case of s H -symmetry A e = —A/, = A. 

If we consider a fictitious homogeneous vector poten- 
tial A, the only correction that we need to make is to 
substitute: 

£n(&) £n (k - ^A^ m ^ n (fe)-%„ fe -A+^2 A ™nk A 

(12) 

Here, we introduced: 



for the group velocity of Bloch electrons, and: 



dki dkk 



(fe) 



(13) 



(14) 



for the components of the inverse band mass tensor. 

The current density operator is obtained by differenti- 
ating Eq. (fTT]) with respect to A. Since we are interested 
in linear response of the current to the vector potential, it 
is sufficient to expand the Hamiltonian up to second or- 
der in powers of A, as we did in Eq. (|12[) . The expansion 
is straightforward, resulting in: 



J= c 



5H 

Ja 



kaB 



E 1 



(15) 



where we introduced the following matrices: 



2 - -i 



(Jk)afS = e (V k )aB ~ — (M k )a[j ' A, 



(16) 



{Vk) a s = <W diag («/,&, -v h -k, v e k, -v e -h) > (17) 

, -H- - o-l o-l o-l \ 

(Mk )*b = 6 a p diag [rn hk , - m h _ fe , m ek , - m e _ fe J . 

(18) 

The average of the single-particle operator Eq. (fT5")) is 
expressible solely in terms of the single- particle Green 



function. However, in the mean-field approximation, the 
corrections to the Green's function itself are expressible 
in terms of the Green's function for the system described 
by Eq. (HH): 



(iu m i - U k ) G k (u m ) = 1- 



(19) 



Expanding (J) to linear order in powers of A, and iden- 
tifying the coefficient of proportionality with the Q re- 
sponse function: 



(J)i = -Q ik A k , 
we arrive at the following expression: 

2 

= ^ E e 



(20) 



Tr ( VifcGfc(w m )VfcfcGfc(u; m ) 



n . h 



(21) 



At this point, we assume a general band dispersion rela- 
tion, with only inversion symmetry £„(— fc) = £ n (fe). Ac- 
cording to Eq. (fl~3|) . and Eq. (fl"4|). we have v n -k = —v n k, 



and m 



7i — k 



b nk- 



The inverse matrix Eq. ([T 



with 



T-Lk given by Eq. (|11[) has a denominator of the form 
(x 2 - E\(k))(x 2 - E 2 _(k))\ x=lUm , where the excitation 
spectrum is given by: 



E±(k) = v (J\m\ 2 + e k ±5 h )* + |A|*. 



(22) 



After taking the traces in Eq. ([21]) , we are left with ratio- 
nal algebraic expressions of x 2 , and the partial fractional 
decomposition is made with respect to this variable. The 
sums of Matsubara frequencies are evaluated by using the 
following identities: 



1 



E 



JbJ m + 



n F (E) - ± 



P ^ {iu m y - E 2 



E 



9(E), (23) 



1 ^ e iUm ° _ g'(E) 

~Ph((^m) 2 ~E 2 ) 2 ~^E-' 



(24) 



where np(E) = (exp(/3i£) + 1) 1 is the Fermi-Dirac dis- 
tribution function. 
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Then, we combine the diamagnetic terms (containing 
factors of the inverse mass tensor) with the "convective 
terms" (containing products of two group velocities) by 
integrating by parts, through the use of the following 
identity which follows straightforwardly from the defini- 
tions Eq. (|T3j). and Eq. (fl4| : 



/ -1\ _1 



dv ink , dv kr , 



dk k 



(25) 



After a simplification, we get the following expression: 



g(E+)-g(E_) ^ , 
± 



E\-E 2 



where, the coefficients L^, and Kik are given by: 

4 Iml 2 A| 2 



(26) 



I 



i k 



\m\ 2 + £ 2 



{Vei ~ V hi )(v e k ~ Vhk), (27a) 



VhiVhk (VH 2 + £ 2 T C) 



(|m| 2 +£ 2 )£± 

+«eiWefc (Vl m | 2 + £ 2 ± + (VhiVek + V ei Vhk) \m\ 



(27b) 



2 2 

e 2 wi 



Equations (|26j) and ([27|) can be simplified further if we 
assume a small deviation from the perfect nesting and 
isotropic bands, defined by Eq. ([3]). In this approxima- 
tion, it is sufficient to assume v e k = —Vhk = Vk-, just as 
in the case of perfect nesting and isotropic bands, which 
give a velocity in the radial direction Vk = Vpk, k = 
{cos(</>), sin(^)}. Here, <fi is measured from the direc- 
tion specified by the nesting vector Q. It may be shown 
that taking into account any deviation in the veloci- 
ties is a small quantity of higher order. Under these 
approximations = (l6|m| 2 | A| 2 /(|m| 2 + £ 2 )) ViV k , 

and {K±) ik = (4|A| 2 £ 2 /(£ 2 + H 2 )£±) W Next, in 
general, the symmetric two-fold 2D Cartesian tensor 
is reducible to an isotropic tensor (proportional to the 
Kronecker symbol Si k ), and a symmetric traceless ten- 
sor. However, one can show that under the above as- 
sumptions, the off-diagonal elements vanish identically. 
Indeed, a reflection relative to the y-axis, for exam- 
ple, changes the azimuthal angle as <f> — > 7r — cf>, and 
cos(2</>) — > cos(2</>), thus the detuning parameter Eq. © 
remains unchanged. However, x — > — x, and y —> y, so 
{xyf{8^)) -> —{xyf(6j^)). Finally, we go from integration 
over k to integration over £, and take the slowly varying 
velocity and density of states in front of the integral as 
their values at FS. We present both the isotropic Q ISO , 
and the anisotropic deviation Q an (so that the principal 
values are Q iso ± Q an ): 



cos(20) 




4|m| 



g(E+)-g(E-) 
E 2 -E 2 _ 



(28) 



Finally, since we are mainly concerned with the pair- 
breaking effect the SDW order has on the Cooper pairs, 
as we vary the detuning parameters in Eq. (|3]), we 
take the T — > limit in the above expressions. From 
Eq. flU]), E q- and Eq. flM]), it follows that g{E) ->■ 

-l/(2E),g'(E) -> 1/(2E 2 ) at zero temperature. 

In order to calculate Q, one needs the values of the 
SC and SDW order parameters |A|, and |m|. These are 
given by the self-consistency equations, that were derived 
clsewhereii for finite temperatures: 



In 



T 



T, 



c.O 



2ttT Rci 



bj m >0 



(E m + i5 k )/E m J_ 
^/{E m +iS^) 2 + |m| 2 u m 



(29a) 



The above expressions arc not defined at zero temper- 
ature. To reconcile with this, we add and subtract 1/E m 
in each of the angle brackets. Then, the combination: 



2nT £ 



uj m >0 



1 

E n 



1 



= In 



ne—> E T 



,T->0. 



We see that the logarithm of temperature cancels on both 
sides of Eq. ([29|) . and the remaining expressions have a 
well-defined zero temperature limit. Using this, and the 
BCS relation between the zero-temperature gap parame- 
ter and the critical SC temperature (|Ao| = 7re _7E T Cj o), 
we get the zero-temperature self-consistency conditions, 
presented above by Eqs. (jj). 



In = 2nT J2 Rc 



1 



T s ,o 



ui m >0 



%) 2 



\ \/{E m t ov^J 

(29b) 

A| 2 . We remind that T C: o is the 
transition temperature for the pure SC state, and T Sj o 
is the transition temperature for the pure SDW state at 



where E m = ^fu 2 



IV. CONCLUSIONS 

To conclude, there is a clear reason behind the sen- 
sitive dependence of the magnetic penetration depth on 
doping in the coexistence phase of SDW and SC orders 
in iron-based pnictide superconductors. In the presented 
model, doping is characterized by a detuning parame- 
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tcr of an electron and hole FSs from perfect nesting, 
which tends to suppress the SDW order in favor of SC 
order, eventually leading to a pure SC state. Although 
the FSs remain gapped due to both order parameters, 
only a fraction of quasiparticles near the FSs contribute 
to the superfluid density, while other quasiparticles par- 
ticipate in formation of the SDW. Therefore, the zero- 
temperature superfluid density increases, as doping is 
increased on the undcrdoped side. This behavior was 
recently observed in experiments! 16 ' 17 Similarly, a pre- 
vious attempt^ to fit experimental data on magnetic 
penetration depth temperature dependence in electron- 
doped BaFe2As2 required an adjustment for the zero- 
temperature value Ao of the magnetic penetration depth 
by factor of 4 between the optimally doped and under- 
doped samples. This adjustment in choice of Ao is con- 
sistent with the suppression of the superfluid density in 
the coexistence region. 

One further prediction that our calculation gives is an 



anisotropic in-plane penetration depth, induced by the 
present SDW order. This would affect the measured 
value for the London penetration depth depending on 
the relative orientation of the sample boundary and the 
SDW nesting vector in mono-crystal samples. In multi- 
crystal samples, the anisotropy is likely to average out. 

We note that although the presented model captures 
the qualitative behavior, there is a need to consistently 
treat the effects of detuning^ and disorder— induced by 
doping not only on the competing orders, but also di- 
rectly on the superfluid density. 
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